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1 .  INTRODUCTION 

The  radar  scattering  propertlaa  of  rectangular  flat  plates  and  rectangular 
cylinders  are  of  concern  for  a  number  of  reasons  of  which  the  following  are 
perhaps  the  most  Important: 

(I)  Some  targets  of  Interest  can  be  approximated  by  one  of  these  shapes. 
In  which  case  the  results  apply  directly, 

(II)  Many  targets  can  be  approximated  by  a  composition  of  these  primitive 
shapes;  calculation  of  the  radar  scattering  properties  of  complex  targets 
My  then  reduce  to  an  iterative  procedure  employing  the  scattering  prop¬ 
erties  of  the  primitives, 

(III)  Since  these  shapes  are  easily  fabricated,  they  are  convenient  for 
anecholc  chamber  experiments;  such  measurements  My  be  used  to  validate 
theoretical  predictions  or  to  assist  In  calibrating  the  anecholc  chamber 
Itself. 

In  this  report,  the  radar  cross  sections  of  flat  plates  and  rectangular 
cylinders  are  calculated  usio&  Physical  Optics,  the  GeoMtrlcal  Theory  of 
Diffraction  and  the  Uniform  Geometrical  Theory  of  Diffraction.  The  theoret¬ 
ical  predictions  are  then  compared  with  experlMntal  measurements  obtained  in 
the  ERL  anecholc  chamber. 

In  the  present  context  the  prlMry  motivation  was  one  of  establishing  the 
domains  of  validity  of  three  important  theoretical  techniques  for  calculating 
radar  cross  sections:  physical  optics,  the  geoMtrlcal  theory  of  diffraction 
and  the  uniform  geoMtrlcal  theory  of  diffraction.  This  objective  Is  part  of 
a  broadar  Investigation  which  siss  at  deriving  guidelines  for  employing  the 
various  technlquss  available  to  the  electroMgnetlcs  researcher.  At  present 
It  appears  to  be  the  case  that  test  methods  yield  sensible  results  over  e 
wider  range  of  situations  Chan  night  be  expected.  In  other  circumstances  the 
opposite  holds.  It  Is  difficult  to  predict  this  behaviour  from  tha  isolated 
theory;  one  must  apply  the  techniques  to  real  probleM  and  confirm  their 
validity  by  comparison  vlth  experlMnt . 

Obviously  even  this  methodology  has  Its  limitations  but  there  Is  no  clear 
alternative  if  one  vlshes  to  be  able  to  select  and  apply  theoretical  Mthoda 
with  confidence  that  the  results  will  ba  Maningful. 

The  rcMlnder  of  this  paper  la  organised  ea  follows.  First  the  method  of 
physical  optlcc  (PO)  is  described  and  applied  in  turn  to  the  flat  plate  ar.d 
tha  rectangular  cylinder.  Next  the  geoMtrlcal  theory  of  diffraction  (GTD) 
and  the  uniform  geometrical  theory  of  diffraction  (UTD)  art  applied  to  the 
same  targets.  A  brief  description  of  the  experlMntal  Masurements  carried 
out  In  tha  anecholc  chamber  la  than  given,  followed  by  a  discussion  of  the 
comparisons  between  theoriaa  and  experlMnt.  Som  limited  comparisons  with 
solutions  obtained  using  the  Mthod  of  moMnta  era  included. 

Finally,  som.'  conclusions  axe  drawn  regarding  tha  dOMlns  of  validity  of  the 
theoretical  techniques. 

2.  SCATTERING  ANALYSIS  USING  PHYSICAL  OPTICS 

2.1  The  Mthod  of  physlcel  optica 

Physical  optica  determines  the  field  scattered  from  a  eeetterer  by  assuming 
that  the  field  on  the  eurface  of  the  scatterer  is  the  geometric  optics 
surface  field.  At  each  point  on  the  geometrically  illuminated  side  of  the 
body,  the  Induced  surface  current  la  tha  same  aa  that  os  an  Infinite 
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tangent  plan*  at  th*  point!  whereas  ovar  the  shadowed  body,  the  surface 
field  la  aero.  For  a  perfectly  conducting  body,  the  surface  current 
distribution  on  the  scatterer  is  then 


J 

■zpo 


2n  *  H1 


in  the  lllualnsted  region, 


In  the  shadowed  raglon 


(1) 


where  n  is  a  unit  vector  normal  drawn  outwarda  f row  the  body  and  g  la  the 
Incident  magnetic  field.  This  la  only  valid  for  wavelengtha  very  much 
smaller  than  the  dimensions  of  the  object.  When  this  surface  current  is 
substituted  into  th*  magnetic  field  Integral  equation,  the  phyalcal  optics 

approximation  for  the  scattered  magnetic  field  Ha  Is  obtained,  that  Is, 


2(n  x  H1)  x  7 


dS . 


(2) 


where  k  -  — ,  R  Is  the  distance  between  th*  field  point  and  the  Integration 
point  on  th*  scatterer  and  I  denotee  th*  Integral  over  the  Illuminated 

S1 

surface  of  th*  body.  The  time  factor  *luC  Is  assumed  and  suppressed. 


When  the  field  point  Is  a  large  distance  from  th*  body,  th*  scattered 
magnetic  field  simplifies  to 


lk 

4s 


-Ikr 


2(n  x  H; 


x  ?  #lht.r’ 


(3) 


where  r  Is  th*  unit  vector  from  th*  origin  to  the  field  point,  r  Is  the 
distance  from  the  origin  to  the  field  point  end  £*  fs  the  radius  vectoj 
from  the  origin  to  the  Integration  point.  For  backseat  ter  ing  r  -  -k  (k 
being  th*  unit  wave  vector),  so  the  backacattered  magnetic  field,  valid  In 
the  far  flald  limit,  la 


—ikr  f  A 

H*  “  — —  j  2(5  x  H1)  x  (-k)  e_ikk'~  dSt  (4) 

JSi 

Useful  accounts  of  the  application  of  phyalcal  optics  to  rudar  acattarlng 
can  be  found  In  reference*  1  and  2 . 

2.2  Sackscatearlng  from  a  rectangular  flat  plata 

k  rectangular  flat  plate  of  side*  2a  and  2b  llaa  In  tba  y-z  plana  aa  shown 
In  figura  1. 
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A  plan*  wave  la  Incident  an  eh*  plat*  In  the  azimuthal  plan*  making  an 

A  A  A 

angle  $  with  the  x-azia.  Thu*  k  •  -eos<>i  -  aln^l  .  Let  the  Incident 

i.  A  ^  / 

magnetic  field  be  H  •  1H  *  .  If  the  incident  wav*  la  a  horizontally 

A  A  a  A  A  A 

pol*ri»«d  vavt,  l  •  n  *  i^,  r*  •  yi^  ♦  *1^  bo  that  this  backac*tt*r«d 
flald  is 


t  Ik  e~lkf 
z  4*  r 


coa*  f  ,21k*^y  dSl 


-  1  £  2H  4ab  coa* 

i  4w  r  o  v  2kasln* 


Hence  the  monoatatlc  radar  croaa  aectlon  of  the  flat  plate  la 


-  C!ltE=iS|2i-i]a  <» 

If  the  Incident  wav*  la  a  vertically  polarlaed  wave,  1  »  ain$i  -  cos$i 
and  the  backscattered  field  la  determined  by  equation  (4)  to  be 


f  ,  11  « 


E  -  (elnd  coadijj  -  cosa4>i  >  ^ 


—  2Ho 


sln(2kasine) 


yielding  a  monoatatlc  radar  croaa  aectlon  equal  to  that  given  by 
equation  (6)  for  horizontal  polarisation.  Although  the  physical  optics 
radar  cross  section  is  polarisation  Independent  In  this  case,  this  is  not 
necessarily  true  of  composite  targets  vhere  multiple  reflections  can  occur; 
this  has  been  strikingly  demonstrated  in  an  analysis  of  the  orthogonal 
trihedral  corner  refl*ctor(r*f , 3) , 

2.3  Backseat taring  from  a  rectangular  cylinder 

A  cylinder  of  height  2b  has  a  rectangular  cross  aactlon  of  aides  2a,  2c. 
The  axis  of  tha  cylinder  lies  along  the  z-axla  as  shown  In  figure  2. 

A  plane  wave  Is  Incident  on  the  cylinder  lt^  the  azijeuthal  plane  making  an 
angle  of  $  with  the  x-axla  so  k  •  -cos^i  -  sin$i  .  Let  the  Incident 

<  A  _4fct  r'  X  y 

■agnatic  tleld  be  Hl  -  1Hq*  ~ 


If  the  Incident  wave  la  horizontally 


polarised,  l  -  1^  and  tha  scattered  field  la  given  by 
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h‘  •  l  f  s -  2H  4b  I*  cos*  «  2ikcco**  ilgl2l£2ls£t 

—  *  4w  r  o  L  2k««ln$ 

,  t  21k«iln4  ■ln(2kcco«4)  I 

+  c  *la*  •  1  ikcCOS*  J 

The  monostatic  radar  cross  asceion  is  then 


0  .  2ikccoe*  aln(2kaalp*) 

°po  Xa  ^  2kaain* 


.  ,2ikasin*  sin(2kccos*>  * 

+  c  Bln*  *  ""itecorT" 


(9) 


This  expression  Is  also  valid  for  a  vertically  polarised  incident  wave  as 
veil. 


3.  SCATTERING  AHALYSIS  USING  KELLER'S  GEOMETRICAL  THEORY  OF  DIFFRACTION 

3.1  The  geometrical  theory  of  diffraction 

In  geometrical  optics  only  incident,  reflected  and  refracted  rays  are 
described.  Keller's  geoaetrlcal  theory  of  diffraction  Is  an  extension  of 
geometrical  optics  in  which  diffracted  rays  are  introduced  by  a  generalisa¬ 
tion  of  Fermat's  principle.  Diffracted  rays  are  produced  when  incident 
rays  hit  edges,  corners  or  vertices  of  scattering  surfaces,  or  when  inci¬ 
dent  rays  impinge  tangentially  on  smoothly  curved  boundaries.  Diffracted 
rays  can  penetrate  into  the  shadow  regions  and  account  for  the  fields 
there.  The  initial  value  of  the  field  on  a  diffracted  ray  la  obtainad  by 
wltiplying  the  field  cn  the  incident  ray  at  the  point  of  diffraction  by  an 
appropriate  diffraction  coefficient.  Diffraction  is  treated  aa  a  local 
phanomanon  determined  entirely  by  the  local  propartiaa  of  the  field,  the 
madia  and  the  boundary  in  the  immediate  neighbourhood  of  the  point  of 
diffraction.  Thua  the  diffraction  coefficient  is  determined  from  the 
solution  of  the  simplest  boundary  value  problems  having  these  local  prop¬ 
erties;  these  problems  are  called  canonical  problems.  Away  from  the 
diffracting  surfaces,  the  diffracted  rays  behave  just  like  tbu  ordinary 

rays  of  geometrical  optica.  Aa  befora,  a  time  factor  i  is  assumed  and 
is  suppressed  throughout  the  work. 

3.2  Rackscattering  from  a  rectangular  flat  plate 

A  rectangular  flat  p'ste  of  tides  2a,  2b  becomes  an  infinite  strip  of  width 
2a  whan  b+“.  Ross(ref,4)  has  shown  that  the  results  for  an  infinite  strip 
can  be  adapted  to  that  for  the  rectangular  flat  plate.  A  place  wave  is 
incident  on  the  flat  plate  in  the  azimuthal  plane  making  an  angle  of  *  with 
the  x-axls  ae  shown  in  figures  1  and  3.  It  can  be  seen  that,  if  corner 
diffraction  is  neglected,  only  edge  diffracted  fields  from  edges  AD  and  BC 
will  contribute  to  the  bsckscattered  field  for  angles  of  incidence  *1*0  as 
for  the  infinite  strip. 

The  bsckscattered  far  field  is  derived  for  an  infinite  atrip  and  the  radar 
cross  section  par  unit  length,  o(length),  la  than  calculated.  Tha  area 
backacattsr  radar  cross  section,  o(area),  for  the  rectangular  plats  ia 
obtained  ueing  the  formula  which  relates  area  and  length  cross 
sections (ref .5) 


-  5  - 
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c(ar*a)  -  -r-  o(laugth)  (10) 

Both  vertical  a  ad  horizontal  polarisation*  are  considered  her*. 

3.2.1  Vertical  polarisation 

Following  Ja**s(r*f.6) ,  the  diffracted  field  for  normal  incidence  of  a 
plane  vav*  on  a  straight  edge  of  a  half-plane  Is 

Ev(sd)  -  D*E*(Pfl)  (ad>  *  e‘lk-*d  (11) 

where  n*  is  the  half-plane  diffraction  coefficient  for  vertical 

polarisation.  E*(Pq)  is  the  incident  field  at  the  point  of  diffraction 

P_  and  *d  1*  the  distance  of  the  field  point  fro*  the  diffraction 

point.  The  diffracted  field  E^  and  incident  field  Ey  have  only 

z-coaponents .  These  diffracted  rays  lie  In  a  plana  perpendicular  to  tha 
edge.  When  a  diffracted  ray  hlta  another  edge  It  can  again  produce 
diffracted  rays  and  so  on.  Thus  on*  has  not  only  singly  diffracted  rays 
froa  tha  two  parallel  edges  contributing  Co  the  backacsttarsd  field  but 
also  doubly,  triply  sod  other  higher  order  sultiply  diffracted  rays 
contributing  **  veil.  The  diffracted  far  field  la  given  by 

equation  (11)  with 

D*(<t>p.$0)  -  ~  l*ac  !<$?  -  pQ)  -  see  J(gp  +  g0>)  (12) 

whera  $o  1*  tba  incident  angle  Measured  fron  the  half  plana  faces  and  <pj. 

la  the  angle  of  the  diffracted  ray  Measured  fro*  th*  ease  half-plan* 
face  at  th*  point  of  diffraction.  Th*  singly  diffracted  ray  fro* 
diffraction  point  A'  la  given  by  equations  (11)  and  (12)  with 

ad  «  p  +  a*  log  where  p  is  the  distance  between  0,  the  coordinate  origin, 
on  the  plet*  and  the  field  point  and  p  »  aaing.  Referring  th*  phase  of 
the  incident  field  to  0,  the  incident  wave  at  A'  la 

E*(A')  -  Eq  e  Th#  air-jly  diffracted  field  scattered  back 

towards  the  source  is  thus 


'  -(»“  I (♦g.  "  +0>  ~  «•«=  *(*A.  *  ♦«)!  E0  faz  (13) 

Figure  3  shows  "  f  “  $•  «o 

<  '  -(*  -  asf) E.  •'!lk**w  ini 

Th*  singly  diffrsetad  ray  fro*  B’  is  given  by  aquation*  (11)  and  (12) 
with  *d  -  p-asing ,  E*(B’)  •  Ep  #ika*la*l  .  I  +  g;  thus 
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*  1  \  .  2ikaeli*  i~ikp 

Ev  '  “v1  +  n*)  Ec  *  .*rn? 


Th«  total  singly  diffracted  field  fro*  A*  and  B’  backacattered  tovtrda 
the  source  la 


cos(2kaalnd)  +  21ka 


eln(2kaelf 


L3Eo 


The  singly  diffracted  ray  fro*  A’  In  the  direction  of  g  ,  •  0,  2n  hire 

d  ^ 

B’  when  a  •  2a  on  the  front  and  rear  of  the  atrip  respectively; 

equations  (11)  and  (12)  constrain  this  field  to  be  aero.  In  order  to 
calculate  the  doubly  diffracted  rays  fro*  B' ,  a  higher  order  edge 
diffraction  ter*  ha a  to  be  considered  -  this  la  called  the  slope  diff- 
reetlon  ter*  where  the  diffracted  field  la  not  proportional  to  the 
incident  field  but  to  the  nor*el  derivative  of  the  incident  field  et  the 
point  of  dlffrectlofi.  This  la  given  by  Ja*ea(ref.6)  to  be 


Ev  *  m  /nkfsla  *<♦? "  *<♦»  -  +„> 


+  ein  1(*?  +  «o)sec*  |(pp  +  (aV*  e‘lk*d  (17) 

for  a  vertically  polstlaed  Incident  wave  at  norval  Incidence  on  the  edge 
of  e  half-plane.  In  calculating  the  Multiply  diffracted  field*,  the 

factor  j  for  grailng  incidence (ref .6)  Is  dropped.  Thla  la  equivalent  to 

adding  the  backacattered  field  resulting  fro*  the  diffraction  of  the 
surface  wave  on  the  front  of  the  strip  to  the  backacattered  field 
resulting  fro*  diffraction  of  the  surface  wave  on  the  rear  of  the  strip, 
both  front  and  rear  waves  contributing  squally.  Thus  the  doubly 
diffracted  field  fro*  B'  In  the  backacetter  direction  is  given  by 
equation  (17)  with  the  normal  derivative  of  the  Incident  field 


Ta  Eo 


-Ikaslng 


+  sit* 


S’  "  2  +  ♦*  *o 


0,  a  -  p  -  aslnd 


whence 


-U2a  -ikp 


ike  /81rk2a 


Similarly  the  doubly  diffracted  field  fro*  A'  le  given  by  equation  (17) 
with  the  normal  darivatlva  of  tha  Incident  field 
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Thus  the  tottl  diffracted  field  is 


_T  .  o  ! 

Ev  "  4a4k4  V 


Eo  (  e  4^aN\  c  j  (1  -  sin4>)  -2iV,a8in<£ 


Cl  4-  8inifr) 
(l  -  sin$) 


Abbs)  /te|j 

wt)  e2ikaein$j 

n<$) 4  *  J 


Higher  order  multiply  diffracted  fields  can  be  obtained  in  a  similar  way 
by  applying  equation  (17)  repeatedly  for  the  various  incident  fields  on 
the  edges  and  summing  the  results.  It  would  be  tedious  indeed  to  have 
to  evaluate  each  higher  order  field  this  way.  It  transpires ,  however, 
that  the  contribution  from  all  the  multiply  diffracted  fields  can  be 
accounted  for  in  a  very  simple  vay  if  one  observes  that  all  multiply 
diffracted  fields  from  doubly  diffracted  fields  onwards  can  be  divided 
into  two  general  cases. 

(a)  The  incident  ray  hits  one  edge  and  then  undergoes  2,  4,  6,  8, 
. . .  diffractions  before  being  backscattered  from  the  other  edge  to 
the  one  it  was  incident  on. 

(b)  The  incident  ray  hits  one  edge  and  then  undergoes  3,  5,  7.  9, 
...  diffractions  before  being  backscatteied  from  the  same  edge  it  was 
Incident  on. 

These  diffracting  mechanisms  are  illustrated  in  figures  4  and  5.  In 
summing  the  contributions  from  these  multiply  diffracted  rays,  it  is 
seen  that  the  doubly  diffracted  field  given  by  equation  (22)  is  a  common 
multiple  in  case  (a)  60  that  the  total  field  fro.-  the  contributions  is 


2E  -ik2a  -ikp 

v  ^  Tka  /8iitk2a  /BlTkp  +  **  +  x  + 


f  1 

ik  2a  2/81 


2v'8frk  ^ 


sin  i(4>  -  )  sec2  J  (♦„  -  *  ) 


rP  e 


P  o' 


d  -iksd 


+  sin  i($p  +  4>o)sec2  l($p  +  i>Q)  j  (sa)  e 


evaluated  at  s  -  2a  and  $p  -  ■»  0,  that  is. 


t  -ik2u 
x  “  4ika  7Biirk2a 


256iirk*a3 


-  9  - 
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Thus  the  sum  Is  the  sun  of  sn  Infinite  geometric  series  and 

-HI  *  Ifo  _i_  (.  + 

v  iks  /8ii?k2a  /8iirkp  cos$  V  256iirk*a9/ 


(31) 


The  triply  diffracted  field  given  by  equation  (27)  Is  a  common  multiple 
in  case  (b)  so  that  the  total  field  can  be  written  as 


E*2 


( 


-ik2a 


-ikp 


(1  -  alnd)  -Zikasin$ 
_(1  +  sin$)a  e 


X  ~  4a^k2  \  /8lTrk2a'y  /81rkp 

(1  ± ,s.lny)  2ikasin<TI  (1  +  *2  +  x-  +  ) 

a  -  sin*)*  e  J  U  +  x  +  x  - ) 


(32) 


where  xz  is  as  given  by  equation  (30). 
Thus 


f2 


4aak2 


(  e~ilc2a\2  e“lkp  f (1  -  sin*)  -2ikasin<|> 

UiSEa)  “ 


(1  +  sin$) 2 


(1.  +  ainj>)  2ifcasin$  ..iC.1.*8  \ 

(1  -  sin*)1  e  J  V  iSbiiTk**1  / 

The  total  backscattered  field  for  the  plate  is  then 


„tot 


^  -  2E 


-ikp 


O  /8iTTkp  Kl 


(33) 


(34(a)) 


where 


Ki 


cos(2kasin0)  +  2ika 


sin(2kasln4i) 

2kasind> 


-ik4a  \  r  -ik2a  . 

e  \  a  1 

256iirk*asy  ltacosi 


-  SPk5 


e  llc2aS\  P (1  -  alnft)  -2ikasin4> 

,'IPvnVia)  jjl  +  sin#)2 


(1+  sln$)  2ikasin$ 
(1  -  8in$)*  e 


(34(b)) 


The  monostatic  radar  cross  section  per  unit  length  for  vertical  polaris¬ 
ation  is  defined  to  be 


ERL-Q344-T* 
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0y(  length) 


11a 


2irp 


,tot 


J*il! 

k 


(35) 


and  the  area  backacatter  radar  cross  section  of  the  rectangular  plate  Is 
given  by  equation  (10)  to  be 


GTD 


¥  Ik, I' 


(36) 


3.2.2  Horizontal  polarisation 

The  diffracted  field  for  normal  Incidence  of  a  plane  wave  on  a  straight 
edge  is 


)  -  D" 


Ej(Po)(sd) 


-i 


-iki. 


(37) 


where  D  is  the  half-plane  diffraction  coefficient  for  horizontal 
polarisation,  E^(Pq)  is  the  incident  field  at  the  point  of  diffraction 
and  sd  is  the  distance  between  the  diffraction  point  and  the  field 
point. 


Dm(<f»p.<f>0)  “  “  l  sec  ”  #„)  +  »ec  H$p  +  $0)]  (38) 


Singly  diffracted  rays  frost  A'  and  B'  are  calculated  in  the  same  way  as 
for  vertical  polarisation  and  the  total  singly  diffracted  field  is 


eh 


cos  (2kasimt) 


8ln(2kasin$) 

21k*  ■  Kail'S*”' 


e-«cp 

v^Irfttp 


(39) 


The  singly  diffracted  field  from  B'  in  the  direction  “  0  hits  A' 
when  s  *  2a.  Thus  the  incident  field  at  A'  that  gives  rise  to  doubly 
diffracted  field  frost  A'  is 


eJ  %  -  2/r  Eo  eik"in* 


-lk2a 


»4?7t3a 


(40) 


The  doubly  diffracted  field  froa  A'  is  given  by  equstlous  (37),  (38)  and 
(40)  to  be 
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4'  '  (4l) 

Again,  the  grazing  Incidence  factor  la  dropped  in  the  derivation  for 

multiply  diffracted  fields  to  allow  for  the  rear  surface  wave  contrlbu- 
t  ions . 

Similarly  the  doubly  diffracted  field  from  5'  can  be  found  to  be 


-v  8  E 


-lk2a  -lkp 


n  ”  o  /8iirk2a  /RwSp 

Accordingly,  the  total  doubly  diffracted  field  is 

.  -ik2a  -lkp  , 

*H  ^  16  So  vfej?  vfep  (43) 

Using  equations  (37)  and  (38)  again  with  the  appropriate  incident  field, 
the  triply  diffracted  field  is  found  to  be 

/  -lk2a\2  -lkp  r  2ikasln4>  -2ikasln$~| 

4  *  *  16  EoV8lTTk2a>  vfep  Li'  -  «ST  *  1  +  sin$  J  <44> 

All  multiply  diffracted  fields  can  be  accounted  for  in  the  same  simple 
way  as  for  vertical  polarisation  (figures  4  and  5).  The  backscattered 
field  for  doubly,  triply  and  all  other  higher  order  diffraction  is  given 
by 


r  -lk2a  -lkp  (  -ik2a\1  -lkp 

s  '  L16  Eo  /fen  /fep  c-3i?  - 16  Eo  V/tan;  .fas  • 

p  2ikasln$  -21kasin$  —] 

|_1  -  sin*  +  1  +  sin*  J(1  +  x*  +  x*  + 


■  [,ec  l(*p  ■  v +  ,ec  +  v1  (,d)  e’lka 


evaluated  at  *  $  *  0  and  a  •  2a.  So 

r  o 


2e~lla» 

/81rk2a 


- 
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This  gives 


„2  _  *_ 


4  *  [ 


e~lk2a  e"1^  1 

^  Eo  /81irk2a  »'BS?p  cos$ 


^  e~ik2ap  #-ikp  |~e2ikasin$  e-2ik*sii4_|  ^  e-ik4a^“ 

EoV/8i7ik2«/  /liirlsp  _  1  -  sin#  +  1  +■  ain$  j  v  41irka ) 


Thus  the  total  backseat tered  field  for  the  plate  for  horizontal  polaris¬ 
ation  is 


stot  a  _  «  « r 

"H  1  —o  kairkp  ‘ 


{49(a)) 


K2  “  coa(2kasin$)  -  2ika 


sla(2kaoln$) 


8  r_i_.  ^p21u*^ , 

/Slirk2a  \  4iuka/  \  cos$  /8itrk2a  j  1  -  sin$  1  +  sin$  I J 


(49(b)) 


The  area  backscatter  radar  cross  section  of  tbe  rectangular  plate  for 
horizontal  polarisation  is  then 


3.3  Backacattering  fro*  a  rectangular  cylinder 

The  rectangular  cylinder  is  shown  in  figure  2.  A  view  of  it  in  tbs 
aziauthal  plane  is  shown  in  figure  6.  Three  diffracting  edges  AZ,  BF  and 
CG  contribute  to  the  backscattered  field  for  the  rectangjlar  cylinder. 
Fro*  the  results  of  the  flat  plate  it  was  seen  that  for  vertical  polarisa¬ 
tion  the  doubly  diffracted  field  which  results  fro*  the  higher  order  slope 
diffraction  ter*  (of  order  k”a)  is  weak  in  comparison  with  the  singly 

diffracted  field  of  order  k  ^ .  Thus  for  vertical  polarisation  only  the 
singly  diffracted  field  needs  to  be  considered.  For  horizontal 

polarisation,  the  singly  diffracted  field  is  of  order  k~^  while  tbe  doubly 

and  triply  diffracted  fields  are  of  order  k~ 1  and  k~  ^  respectively.  Hence 
multiply  diffracted  fields  need  only  be  considered  for  borisontal  polarisa¬ 
tion.  The  technique  used  to  handle  multiply  diffracted  field*  in 
Section  3.2  proved  rather  cumbersome  and  would  be  even  aore  so  for  tbs 
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rectangular  cylinder.  Rudduck(ref .7)  has  pointed  out  a  *nch  eore  elegant 
way  of  computing  the  nett  multiply  diffracted  field.  This  approach,  known 
as  Che  self-consistent  field  method  is  employed  here  to  calculate  the 
Multiply  diffracted  field  for  the  case  of  horizontal  polarisation. 

3.3.1  Vertical  polarisation 

The  diffracted  field  for  normal  Incidence  of  a  vertically  polarised 
plane  wave  on  the  straight  edge  of  a  vedge  of  angle  ^  ia(rtf .6) , 


Ey(sd)  D*  Ej(Po)(sd)  1  e"ik*d  (51) 

where 

D*(V+o>  ■  !$k  [[co*  T"  cos  §(♦?  -  v] 

-  jjeos  3  -  e«s  3<^p  +  *0>J  j  (52) 

is  Keller's  diffraction  coefficient  for  vertical  polarisation  for  a 
vedge  of  angle  Is  the  incident  angle  Measured  from  one  of  the 

wedge  faces  and  $p  is  the  angle  of  the  diffracted  ray  measured  from  the 

same  wedge  face.  Referring  to  figure  6,  the  singly  diffracted  far  field 
from  A1 ,  B*  and  C'  can  be  calculated  using  equations  (51)  and  (52)  with 
the  following  information  for  A’ ,  B'  and  C’ : 


TABLE  1.  EXPRESSIONS  FOR  INCIDENT  FIELD,  DIFFRACTING  ANGLE  AND  DISTANCE  TO 
THE  FAR-FIELD  POINT  FOR  DIFFRACTING  POINTS  A' ,  B’  AND  C' 


D 

sd  ■  p  +  asin$  -  ccos$ 

♦a*  ■  *o  "  i  •  * 

Ei(A')  -  E  .-iMMi«4-ee»«+) 

0 

B' 

sd  ■  p  -  asin$  -  ccos$ 

V  -  *0  -  f  +  * 

Ei(B')  -  E  eiV(asin^fccos^) 

0 

C’ 

sd  -  p  -  aaln$  +  ccos$ 

♦c'  “  *o  ‘  * 

Ei(C')  -  Eq  eik(asinP-ccos«) 

The  phase  of  the  Incident  wave  is  referred  to  the  coordinate  origin,  0. 
p  is  the  distance  between  0  and  the  field  point  and  for  far  field, 
p  »  aslnd,  p  »  ccos$.  The  total  singly  diffracted  field  for  vertical 
polarisation  is 

r  .  M  p  F 
^v  3  o 


(53) 
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,  ,  2ikccos$  -2ikasind 

_  4  2ikccos*  .  e  T  e  T 

F  -  _  _  cos(2kasin$)  e  T  +  - r-— - - - 

2sin  sin  ^(ir  -  $) 


e2ikasin$ 
in  |(tt  +  $) 


tJ  21k(asln$-ccos$)  2 

T  C  —  T  ' 


2sin  i(*  +  2$)sln  |(ir  -  2$) 


the  monostatlc  radar  cross  section  per  unit  length  for  vertical  polaris¬ 
ation  Is  given  by 


E  a 

Ov(length)  -  2rp 


so  the  area  aonostatlc  radar  cross  section  is 


ycn>  .  |bi  lF|l 
v  3r  ■  • 


3.3.2  Horizontal  polarisation 

The  diffracted  field  for  uoraal  incidence  of  a  horizontally  polarised 
plane  wave  on  the  straight  edge  of  a  wedge  of  angle  j  is (ref. 6), 


Eg(sd)  -V  D*  eJ(Po)  (Sd)  J  e~lk*d 


\  2/T  i  ;  2tt  1.  ,  T)-1 

D  (VV  "  sTST^'X008  “  "  c  •  3(*p  •  *o^j 

j”“  — )  ■—  l"' 

+  jcos  ^  -  cos  +  i0)J  J  (58) 

is  Keller's  diffraction  coefficient  for  horizontal  polarisation  for  a 
wedge  of  angle  j.  The  total  backscattered  field  can  be  represented  as 

the  sun  of  the  singly  diffracted  field  and  the  multiply  diffracted 
fields. 

(a)  Singly  diffracted  field 

Referring  to  figure  6  the  singly  diffracted  field  Is  calculated  from 
equations  (57)  and  (58)  for  k' ,  B' ,  C*  with  Table  1  giving  the 
required  Intonation.  Thus  the  total  singly  diffracted  field  for 
horizontal  polarisation  la 
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(59(h)) 


-  |  e21kccos0  COB(2kaain6)  +  £. 


21kccos4>  [  2ikasin$ 


2»in  Nln  |(tt  +  $) 


_  e  2^-lutg:(p4l  ^  #2ik(asind-ccosd)  2  +  _ 1 

sin  2(*  -  $)  5  2n In  jC*  +  2$)  sin  i(it  -  2$) 


lln$ 

-  ♦)_ 


(59(b)) 


(b)  The  aelf-consistent  GTD  formulation  for  the  multiply 
diffracted  fields 

The  multiply  diffracted  fields  result  from  waves  travailing  at  least 
once  along  the  surface  of  the  cylinder.  There  are  infinitely  many  of 
these  waves  on  each  'ace  of  the  cylinder.  The  sell-conslstast  CTO 
concept  combines  all  the  waves  on  each  face  into  two  waves  travelling 
in  opposite  directions  with  unknown  complex  amplitudes  as  shown  in 
figure  7. 

The  complex  amplitudes  Ct,  C2 ,  ...  C,  are  determined  by  a  self- 
consistent  field  procedure.  Eq,  the  amplitude  of  the  incident  wave, 

is  factorised  out  so  that  the  unknown  {C^}  result  from  a  unit 
amplitude  incident  vave. 

At  A’  * 


,  -ik2a  ,  \  -lk2a 

C,  -  Vj  +  i  Cs  D*(0,  0)  *-7fj-+  5  C„  O*  Oj  (60(a)) 

V i  represents  the  direct  source  contribution  to  Ci,  that  is,  the  wave 
diffracted  from  B'  to  A'  for  a  unit  amplitude  wave  incident  on  B1, 
the  second  term  represents  ths  contribution  from  the  wave  travelling 
along  surface  A'  -  B'  being  diffracted  at  B*  to  A'  and  the  third 
term  gives  the  contribution  from  the  vave  travelling  along 
surface  C'  -  B*  being  diffracted  et  B’  to  A'.  The  iiffraetion 

coefficient  la  multiplied  by  for  gracing  incldence(ref  .6) . 
Similarly, 

.  -ik2c  ,  \  -lk2c 

C,  -  V,  +  |  C,  D*(0,  0)  +  }C,  (60(b)) 


At  B* 

1  .  a"lk2#  1  «/3v  \rikU 

C,  *  \  C,  0*(0,  0)  *^=-  ♦  ±  Cs  (60(c)) 
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.  -ik2c  /  \  -Ik2c 

c*  -  V*  +  i  Cs  D*(0,  0)  *7^-  +  ic«  Dm^j,  o J  *7——  (60(d)) 


At  C', 


-lkJtc 


Cs 


.  -ik2c  ,  a  v  -lit 

-  Vs  +  j  C*  D*(0,  0)  *7^-  +  i  C,  dP(*§,  Oj  *7£ 


(60(a)) 


Cs 


1  .  #'lk2*  1  -A*  N  ,~ik2» 

Vs  +  i  C7  Da(0,  0)  *7—-  +  5  C.  Da(^%  oy  *757—  (60(f)) 


At  D', 


.  -Ik2a  . 

C7  -  V7  +  5  Cs  D*(0,  0)  C7—  -  +  J 


\  ~lt2» 

Cs  "  2*  °j  (60(g)) 


C.  -  V.  +  i  Cs  D*(0,  0)  *7^—  +  J  Ci  D*(^~,  o) 


-lk2c 


(60(h)) 


The  waves  excited  by  the  Incident  wave  are  given  by 


V< 


eikCa.in*+ccos*)  ^  *  +  ^  C^l 


vj  -  0 


(6) (a)) 
(61(b)) 


Vi 


-ik(asin6-ccoa6) 


/  \  -lk2a 

D-(0,  2  ~  )  V5T“  (61(c)) 


Vs 


-  eik(«*l^ccna«)  j,a(0>  ^ 


-ik2c 

*i^2c 


(61(d)) 


Vs 

Vs 


s  ik(aaind+ccoed)  _» 


-lk2c 


(°,  w  -  ♦)  *7^— 


(61(e)) 

(61(f)) 


V7  -  .ihUeind-cco^)  ^  3*  _  ^ 


t-lk2a 

r2m 


(61(g)) 


-ik2e 


V.  -  B-lk(aaln#-cc.o#d)  0.(0i  „  +  4) 


(61(h)) 
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The  unknown  complex  amplitudes  can  then  be  solved  bp  matrix 
Inversion, 


C  -  t"!V  (65) 


The  total  backscattered  field-  El!,  from  these  waves  diffracted  at  A' , 
B'  and  C'  la 


<  -  *.  ^ 


(66(a)) 


where 


G2  -  |  j^jcj  d"^|  -  $,  oj  +  Ci  D*(ir  +  $,  0)] 


-ik(asinp-ccoa$) 


+  fCi  „■(*  +  0^  +  Ca  D*(7T  -  ♦.  0)1  elk(Mln^CC0^) 

L  v  /  _j 


+  [  Is  D !*(♦,  0)  +  Ct  D*(j|  _ 


ik(aBind-ccos$) 


J 


(66(b)) 


(c)  The  total  backscattered  field 

The  total  backscattered  field  is  defined  by  E^0t  -  Ejj  +  E^,  that  is. 


mr  «"ikpr  G» 

s  ^  Eo  *7rlsm  + 


C2] 


(67) 


The  area  monostatic  radar  cross  for  horizontal  polarisation  can  then 
be  written 


*T°  -  8b>k 


Gi 

/6lwk 


+  Gj 


(68) 


4.  SCATTERING  ANALYSIS  USING  THE  KOUYOUKJ IAN-PATHAK  UNIFORM  GEOMETRICAL 

THEORY  OF  DIFFRACTION 

4.1  The  uniform  geometrical  theory  of  diffraction 

GTD,  aa  originally  developed  by  Keller,  falls  in  the  transition  regions 
adjacent  to  shadow  and  reflection  boundaries  where  the  solution  becosMs 
infinite  and  discontinuous.  Kouyoumjlan  and  Pathak(ref  .8)  modified 
Keller's  diffraction  coefficients  in  the  GTD  solution  by  s  Multiplication 
factor  involving  a  Fresnel  integral.  Whan  the  observation  point  is  close 
to  a  shadow  boundary,  the  modification  factor  approaches  saro  so  as  to 
compensate  for  the  singularity  in  Keller's  diffraction  coefficients,  and  to 
give  rise  to  a  finite  field  solution  there.  They  refer  to  this  approach  as 
the  uniform  CH>  or  DTD,  «  name  adopted  here  to  distinguish  it  froai  GTD  with 
Keller's  diffraction  coefficients. 
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For  normal  Incidence  on  a  wedge  with  exterior  angle  nr,  the  Kouyounjlan- 

Pathak  diffraction  coafflclanta  D*,  D*  for  vertical  and  horitontal  polar¬ 
isation  respectively  are 

De-a,  v  V  -  -  fj£j£  [cot(!+r)  F(kU+(B*)) 

+  cot(*  ~n-"j  P(kLa‘(B“)>  +  j^cot(i-y^  F(kLa+(B+)) 

+  cot^-frft  J  F(kU-(B+))J  j  (69) 

where 

B‘  -  dp  -  d0. 

a+  - 

dB>  d  as  defined  In  Section  3. 

r  O 

F(x)  -  21/jT  e1*  f  e"lTJ  dt;  (70) 

a±(B)  -  2  cos*  — (71) 

In  which  N*  are  the  integers  which  most  nearly  satisfy  the  equations 

2nn  N+  -  B  »  n, 

2im  -  B  - 

parameter  given  by 

for  plane-wave  normal  Incidence 

for  cylindrical-wave  Incidence  (73) 

for  conical-  and  spherical-wave  normal  Incidence 


l  is  the  distanci 
r 


|  r  +  t' 

d  1 

a  a 


(72(a)) 

(72(h)) 
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where  ei  la  che  dlaceace  of  the  source  point  to  the  diffraction  point,  ad 
la  tha  diatance  of  the  field  point  to  the  diffraction  point,  r'  la  the 
radlue  of  tha  incident  cylindrical  wave  and  r  la  the  perpendicular  diatance 
of  the  field  point  fro*  tha  edge. 

If  tha  arguments  of  tha  four  tranaltlon  functions  in  equation  (69)  exceed 
MO,  the  tranaltlon  functions  are  approximately  equal  to  one  and  equa¬ 
tion  (69)  reduces  to  Keller's  diffraction  coefficients. 


4.2  Backscatterlng  frcm  a  rectangular  flat  plate 

The  UTO  backacattered  radar  cross  section  of  the  rectangular  flat  plate  la 
derived  here.  As  discussed  in  Section  3.3  multiply  diffracted  fields  are 
only  significant  for  horixontal  polarisation.  In  handling  the  multiply 
diffracted  fields  for  horltontal  polarisation,  the  sore  elegant  self- 
consistent  field  method  la  eaploysd. 


For  diffraction  from  a  half-plans,  n  •  2  and  che  diffraction  coefficients 
in  equation  (69)  reduce  to 


D*'"a.  ♦p,  do)  - 


-i*Mp 

27sr|j 


F(kLa(8  ))  r  F(kLa(B 

*  1  X  ■ 


cos(6_/2)  cos 


a(S*))~ 

(P+/2)_ 


(74) 


where 


ai(8)  -  a<8)  -  2  cos* (8/2). 


4.2.1  Vertical  polarisation 

The  multiply  diffracted  fields  are  negligible.  The  singly  diffracted 
fields  from  edges  A'  and  B'  (figure  3)  art  given  by  equation  (11)  with 
D  as  defined  in  equation  (74) .  The  distance  parameter  L  for  plans  wave 
Incidence  and  far-field  backacattered  field  is  very  large  so  that  the 
KouyouaJ lan-Pathak  diffraction  coafficlent  reduces  to  that  of  Keller. 
The  singly  diffracted  field  is  thus  given  by  equation  (16)  and  the  area 
backscatter  radar  cross  section  la 

_UTO  4b*  |  .  ....  aln(2kaein4) 

o  -  -5-1  co*(2kasln4)  +  21ka  — stL-i’.-tr— 

V  » 


(75) 


4.2.2  Horixontal  polarisation 

(a)  Singly  diffracted  field 

Aa  in  Section  4.2.1,  the  singly  diffracted  field  Is  the  same  as  that 
obtained  with  Keller's  diffraction  coefficient.  Thus  the  singly 

diffracted  field  e|  Is  given  by  equation  (39). 

(b)  The  self-consistent  DTD  formulation  for  the  multiply  diffracted 
fields 

The  diffracted  waves  travelling  In  opposite  directions  with  unknown 
complex  amplitudes  C>  C)  Kq  on  one  side  of  the  strip  and  Cj  Eq, 

C*  Eo  on  the  back  side  of  the  strip  ars  illustrated  in  figure  8. 
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At  B' 


/  \  -ik.2*  r  — 

Cl  -  ,»-lkaBini*  Dm(^2a,Of2  -  +  \  C,  Dm(a,0,0)  +  C*  DB(a,2it,0) 


fla 

(76(a)) 


C2  -  e_lkasin*  D*(2a,2irtI  -  ^)*^— +  |  C„  D*(a,0,0)  +  C,  Dm(a,2Tr,0) 


-  e-i^2a 


~7TT~ 

(76(b)) 


At  A’, 


ikasin*  _®/l  .  it  .  Ae  1 

D  ^2a*°*2  +  V^TT+  1 


Ci  •  s 


Cx  D  (a,0,0)  +  C2  D 


i®(a,2it,0)  - 


-ik2a 


72a 

(76(c)) 


C%  “  c 


ikasiniti 


/  v  -iV2a  i-  -|  -J-k- 

Dmf 2a,2rr,|  +  4>j  £^=  ~  +  |  C2  D*(a,0,0)  +  C,  D  (a,2ir,0)  —fg* 


-ik2a 


(76(d)) 


From  equation  (74) ,  it  is  seen  that 


DK(a,0,0)  -  -D*(a,2n,0) , 

Dm(^2a,0,I  -  $)  -  -D®(^2a, 2u , j  -  $)  and 

0Bi2aI0f|  +  ■  -DB^2a.2ir.|  +  . 


Substituting  these  into  equation  (76)  gives  Cs  -  -Clt  Cj  -  -C»  and 
the  following  equations  for  Ci  and  C3 > 


(f  v  —  ik2a  -am 

2a,0,|  -  ^  £7ff—  +  C3  D“(a,0,0,)  s-^ 


-ik2a 


(77(a))  * 


/  -ik2a  -ik2a 

C,  -  D*(2a,0,5  +  i)  S,^T~  +  ci  D®(a,0,0)  6/^" 


Solution  of  equation  (77)  yields 


X  +  ZJ 

Cl  "  ttf 


(77(b)) 


(7R(a)> 


f 
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Y  +  ZX 


1  -  Z1 


(78(b)) 


/  \  -ik2a 

X  -  e-lka8i^  D^a.O.f  -  «)  £^r' 


Y  -  elka8ln*  DBf2a.0.|  +  ^) 


-ik2a 

Z  ■  D°(a,0,0)  • 


(78(c)) 


(78(d)) 


(78(e)) 


Fro*  B' ,  the  multiply  diffracted  field  in  the  backscattered  direction 
is 

/  \  s  -ik(p-a»in$) 


E^1  *  \  rClEo  D°(2a.2  ♦  *.0)+  C8Eo  DB^2a.-2  -  *.o)J  ^ 


(79(a)) 


and  from  A*  > 


-ik(p+asin<|>) 


/  .  /  v~l  -lKVp-t-aa 

eJ2  -v-  ~  [c,Eo  D“(2a.|  -  4>,0j  +  C*Eo  D*(2a,-2|+  $.0jJ  5  /= 


(79(b)) 


But  Ci  »  -Cx,  C,  -  -C»  and  equation  (74)  gives 

0B4a,I  +  *.o)  -  -D"(2a.%  -  <i»»0 


D°^2a,“  -  4,0 j  -  -B*^2a.^'+  <t»0y 


Ejl  -u  CtEo  Dffi^2a,^  +  eita8ltl4> 

*5?  -v  C,5w  E*(2a,I  -  ♦.o)*7F-.  tVa8in* 


(80(a)) 


(80(b)) 


mmss&sm  .  yJt. i- 'J  i.i ; .i  ;i '•«. i;-».  X" *  % 


-  23  - 


ERL-0344-TR 


The  total  backacattersd  field  from  multiply  diffracted  rays  la 


4 


D-  2a,, 


+ 


eikasin$ 


+  Cj  D 


e-ikasin$ 


(81) 


(c)  The  total  backscattered  field 

The  total  backscattered  field,  Elot,  is  the  sum  of  the  singly  and 
aultiply  diffracted  fields, 

-lkp 

^  ^  Eofi^“«  <82<a» 

where 

"  -  - 

♦  C,  D-(2,.I  ,  *.o)  -  C,  -  *.o)  (82(b)) 

Ci  and  Cj  are  given  by  equation  (78) . 

The  area  backacattcr  radar  cross  section  for  horizontal  polarisation 
is  then  glvan  by 


_UTD 


8kbz | Q | * 


(83) 


4.3  Backscattering  from  a  rectangular  cylinder 


The  UTD  backscatter  radar  cross _  sections  for  the  rectangular  cylinder 
(figures  2  and  6)  are  derived  here.  The  derivation  follows  that  of 
Section  3.3  closely  with  the  diffraction  coefficients  of  Kouyoumjian  and 
:-athak  (equation  (69))  used  instead  of  those  of  Keller.  The  wedge  angle  in 

this  case  is  j  so  that  Wf  »  — . 

4.3.1  Vertical  polarisation 


Only  singly  diffracted  field  contributions  are  significant.  As  in 
Section  4.2,  the  Kouyoumj ian-Pathek  diffraction  coefficient  reduces  to 
the  Keller  diffraction  coefficient  for  far-field  backscattered  singly 
diffracted  fields  so  that  the  backscattered  radar  cross  section  is  the 
same  as  that  given  by  equation  (36),  that  is. 


a0TD 

v 


(84) 


with  P  defined  in  equation  (34) 
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4.3.2  Horizontal  polarisation 

The  singly  diffracted  field  is  again  the  same  as  that  derived  using 
Keller's  diffraction  coefficient  and  is  given  by  equation  (59). 

The  multiply  diffracted  fields  (figure  7)  are  as  given  by  equation  (65) 
with  the  elements  of  matrices  V  and  T  now  defined  by 

V,  -  e^Cln^ccoa*)  D-(2a,0,*  +  *)  ,  (85(a)) 


v2  -  0, 


(85(b)) 


-lk(asin$-ccos$) 


(za.0,1  - 


m  * 


(85(c)) 


V,  -  aik(“la*-ccoe«  D*(2c.O.*) 


(85(d)) 


— it2c 

V5  _  elk(asin^cco.^)  D»(2c>0>lr  .  ,  (85(e)) 


v6  -  0. 


V7  -  e 


ik(asinif-ccos$)  m 


(2a, 0,^|  -  <)>) 


(85(f)) 


(85(g)) 


V,  -  e-ih(aHin$-ccos(0  D»(2C(0flr  + 


(85(h)) 


.  -Ik2a 

|  Da(a.  0,  0)  *75- —  , 


(86(a)) 


T2  *  2  D"(rfV  %  °)  fi7?r1 


(86(b)) 


T3  •  ^  D  (c.  °» 


(86(c)) 


T4  -  I  D •(-&-  SL  Q\ 

T4  2  D  ViTTT'  2*  7  /ft  * 


(86(d)) 
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The  multiply  diffracted  field  in  the  backscattered  direction  Is 


-v  E  s~r=~  <5, 
o  /p  s 


(87(a)) 


G,  -  |  [  [c,  D“(2a.|  -  *,o)  +  C2  +  ♦,0>]  e-^(«l^-cco^) 

+  [c,  D*(2.,|  t  *,o)  +  C*  DB(2c,ir  -  *,of|  .ikCa.infrtcco.*) 

+  jj.*  d“(2c,$,0)  +  C*  Dm^2a,^j-  ♦,0^1  elk<“ln*“CCOB*)j  (87(b)) 


The  total  backscattered  field  la  thus 


-tOt  *  * 

hi  '  E 


-ikp  |  G 


o  /p 


~  Gj  “I 

yzm  +  G» 


where  Gt  Is  as  given  In  equation  (59(b)), 

The  area  mono static  radar  cross  section  for  horizontal  polarisation  Is 
then 


5.  EXPERIMENTAL  MEASUREMENTS 
5 . 1  The  anechoic  chamber 

The  experimental  observations  carried  out  to  check  the  theoretical  predic¬ 
tions  were  performed  in  the  DRCS  anechoic  chamber.  This  facility  has  a 
normal  source-target  separation  of  12  a,  s  frequency  range  extending  from 
2,5  GHz  to  40  GHz,  and  analogue  recording  of  neasurenenta(ref .9) , 

A  problem  that  has  to  be  combatted  in  the  anechoic  chamber  measurements  is 
null  drift,  that  is,  a  drift  in  the  adjustment  for  background  cancellation. 
This  can  arise  from  frequency  instability  and  temperature  variation  In 
passive  components.  The  low  power  returns  are  especially  susceptible  to 
any  null  drift  that  occurs  as  Is  evident  In  the  high  incidence  of  data 
aaynaetry  in  the  flat  plate  results,  particularly  at  the  lower  frequencies, 
and  in  the  10  cm  cylinder  at  4  GHz  for  vertical  polarisation. 

In  general  the  performance  of  the  chamber  Improves  with  increasing  fre¬ 
quency  because  the  antenna  gains  are  higher,  the  antenna  beamwldths  are 
narrower  so  that  less  of  the  wall  is  illuminated  and  the  radar  absorbent 
material  is  more  effective.  However,  there  are  some  isolated  narrow 
bandwidths  within  the  operating  frequency  range  for  which  the  chamber 
performance  is  uncharacteristically  poor.  The  12  GHz  frequency  has  been 
Identified  as  one  of  them;  a  dramatic  increaae  In  background  noise  being 
experienced  at  the  time  of  measurements  at  this  frequency. 
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Apart  fron  the  poor  chamber  performance  another  contributing  factor  to  tba 
unreliable  reaulta  at  12  GHz  could  be  the  frequency  instability  of  the 
Klystron  oscillator  used.  Frequency  instability  increases  with  frequency 
and  will  result  in  null  drift  problems .  Sobs  of  the  10  GHz  results  see*  to 
be  similarly  affected,  the  10  GHz  and  12  GHz  meaaurementa  are  done  using 
the  same  Klystron  oscillator  while  the  4  to  8  GHz  Measurements  are  done 
using  another  . 

As  the  operating  frequency  Increases  other  problems  are  decreased  sensitiv¬ 
ity  of  the  receiver  and  components  and  greater  attenuation  in  the  cables. 

5.2  Flat  plate  experiments 

The  flat  plate  measurements  were  conducted  using  a  brass  rectangular  plate 
10  cm  x  10  cm  and  thickness  1  mm  (figure  1).  Five  frequencies  were 

employed,  viz  4,  6,  8,  10  and  12  GHz,  corresponding  to  values  of  as 

shown  In  Table  2.  Measurements  were  made  over  the  range  -90*  to  490*  so 
that  symamtry  could  be  used  to  assess  system  errors.  The  analogue  X-»Y 
plotter  output  was  sampled  at  2*  Intervals  for  analysis. 


TABLE  2.  DIMENSIONS  OF  FLAT  PLATES  IK  TEEMS  OF  WAVELENGTH 


Frequency 

(GHz) 

4 

6 

8 

10 

12 

Ratio  of  aide 
of  plate  to  , 
wavelength,  — 

1.3 

2.0 

2.7 

3.3 

4.0 

5.3  Rectangular  cylinder  measurements 

The  method  adopted  in  the  theoretical  analysis  uzlng  the  Geometrical  Theory 
of  Diffraction  ignores  diffraction  from  the  ends  of  the  cylinder.  In  order 
to  assess  the  significance  of  this  effect,  three  cylinders  of  equal  cross 
sections  10  cm  x  10  cm  but  with  lengths  10,  30  and  50  cm  were  used.  The 
cylinders  were  fabricated  of  brass  of  thlcknsss  1  mm  to  s  nominal  tolerance 
of  0.5  mm. 


Again,  five  frequenciee  were  employed,  namely,  4,  6,  8,  10  and  12  GHz.  The 
dimensions  of  tha  targets  In  tarns  of  Incident  wavelength  are  given  In 
Table  3.  Measurements  were  made  over  the  range  -90*  to  490*;  tha  analogue 
X-Y  plotter  output  was  sampled  at  2*  intervals  for  analysis. 


*.  The  frequency  stability  has  sines  been  Improved  by  phase  locking  tha 
Klystron  to  a  atabla  source. 
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TABLE  3.  DIMENSIONS  OF  RECTANGULAR  CYLINDER  IN  TERMS  OF  WAVELENGTH 


Frequency 

(GHz) 

Ratio  of  width  of 

cylinder  to 

wavelength, 

ilatio  of  length  of  cylinder 
2b 

to  wavelength, 

I.etgth  of  cylinder, 

2b  (cm) 

10 

30 

50 

4 

1.3 

1.3 

4.0 

6.67 

6 

2.0 

2.0 

6.0 

10.0 

8 

2.7 

2.7 

8.0 

13.3 

10 

3.3 

3.3 

10.0 

16.7 

12 

4.0 

4.0 

12.0 

20.0 

S.4  Asymmetry  of  ante hole  chamber  measurements 

One  of  the  features  of  the  anecholc  chamber  measurements  presented  In 
figures  9  to  48  Is  the  occasional  asymmetry  between  the  0  to  90*  data  and 
the  0  to  -90*  data.  The  models  were  not  of  high  quality  machine-shop 
standard  so  a  certain  amount  of  asymmetry  was  anticipated  but  several  of 
the  RCS  patterns  for  the  flat  plate  show  Inexplicable  variations.  Possible 
explanations  for  these  variations  Include  (1)  model  construction  errors , 
such  as  curvature  resulting  in  focussing  effects.  (11)  errors  In  mounting 
the  models  vertically  In  the  chamber,  leading  to  sampling  over  oblique 
arcs,  (111)  errors  In  orienting  the  model  at  $  •  0*  leading  to  a 
translation  of  the  pattern  ocross  $  (see  figure  31),  (lv)  null  drift 
problems  discussed  In  Section  5.1  end  (v)  errors  Introduced  by  the  snalogue 
plotting  equipment. 

Of  these,  (v)  is  bslieved  to  be  most  unlikely  while  the  effect  of  null 
drift  Is  expected  to  be  more  significant  at  the  lower  f requencias ,  4  to 
6  GHz,  when  power  returns  are  low.  (ill)  could  explain  snail  deviations 
but  not  the  larger  ones.  (1)  would  ba  an  appealing  solution  wars  It  not 
for  the  feet  that  for  both  fist  plats  and  cylinder,  the  errors  srs  greatest 
at  4  GHz,  decreasing  through  6  GHz  to  only  <v2  dB  for  most  maaaursments  at 
8  GHz  and  10  GHz,  but  Increasing  again  at  12  GHz.  If  the  12  GHz  results 
are  discounted  because  of  the  poorer  performance  of  the  chamber  at  12  GHz, 
nsither  (1),  (11)  nor  (111)  would  seem  to  account  for  the  asymmetry  because 
errors  in  6  and  9  would  be  much  more  critical  owing  to  the  finer  lobe 
structure  at  higher  frequencies.  Thus  ths  most  likely  explanation  of  the 
discrepancies  is  unidentified  problems  with  the  anecholc  chamber  Itself. 

Unfortunately  other  demands  on  ths  anecholc  chamber  resources  mads  a  second 
aeries  of  measurements  out  of  ths  quastion  so  the  analysis  has  bean  under¬ 
taken  In  the  knowledge  that  the  deviations  from  symmetry  may  wall  reflect  a 
variety  of  errors  as  discussed  above.  The  majority  of  the  measured  data 
shows  no  significant  asymmetry  (compared  with  the  variations  of  RCS  with 
9) ,  however,  so  there  saems  to  be  a  strong  case  for  retaining  the  exper¬ 
imental  points  for  comparison  with  tbs  thaorseicsl  results. 
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6.  RESULTS 


6.1  Flat  plate 

6.1.1  Vertical  polarisation 

At  4  GHz  the  PO,  GTD,  UTD  and  anechoic  chambc-  (AC)  estimates  agree  to 
within  1  dB  out  to  about  IS*  from  normal  incidence,  after  which  the 
theoretical  values  depart  grossly  from  the  measured  RCS  (figure  9).  The 
AC  values  are  symmetric  within  ±2  dB,  increasing  confidence  in  their 
accuracy.  The  6  GHz  AC  measurements  are,  by  contrast,  in  fair  agreement 
with  PO,  GTD  and  UTD  out  to  45*  (figure  10).  Beyond  45*  the  0  to  90*  AC 
values  diverge  significantly  from  the  0  to  -90*  AC  values,  differing  by 
14  dB  at  62*.  The  UTD  estimates  near  grazing  incidence  agree  best  with 
AC  values,  PO  grossly  underestimates  it  while  GTD  is  singular  there. 

The  agreement  between  AC,  GTD  and  UTD  values  at  8  GHz  is  excellent  out 
to  55*  where  measurements  show  a  shallow  minimum  which  is  not  predicted 
by  GTD,  UTD  or  PO  (figure  11).  Again  UTD  succeeds  best  near  grazing 
incidence  while  PO  and  GTD  exhibits  the  same  problems  as  at  6  GHz.  At 
10  GHz  agreement  between  UTD  and  AC  is  within  ±2  dB  over  the  entire 
range  (figure  12)  while  GTD  only  falls  near  grazing  incidence.  The 
12  GHz  results  (figure  13)  are  in  conflict  with  the  trend  shown  from  4 
to  10  GHz.  As  mentioned  in  Section  5.4,  the  12  GHz  AC  data  is  unrelia¬ 
ble  and  so  this  deviation  from  the  observed  trend  nay  be  an  artifact. 

In  all  cases  reviewed  above,  PO  consistently  underestimates  the  RCS  for 
aspect  angles  away  from  broadside,  in  uddition  to  predicting  unphysical 
nulls.  The  GTD  and  UTD  solutions  are  identical  up  to  within  10*  or,  in 
some  cases  5*  of  grazing  incidence.  The  infinite  fields  predicted  by 
GTD  at  grazing  incidence  are  unpbyslcal  whereas  thoss  predicted  by  UTD 
remain  finite  and  agree  with  AC  values  to  within  f3  dB  at  6,  8  and 
10  GHz. 


6.1.2  Horizontal  polarisation 

PO  solutions  for  convex  tsrgets  are  independent  of  polarisation;  this 
is  not  the  care  for  GTD  and  UTD. 

The  4  GHz  results  show  an  agreement  between  theory  and  experiaumt  out  to 
18*,  beyond  that  to  70*  the  GTD  and  UTD  solutions  mstch  the  AC  values  tc 
within  *8  dB  (figure  14),  Beyond  70*  UTD  and  PO  values  decrease  very 
rapidly  while  AC  values  remain  fairly  constant  and  GTD  becomes  singular 
at  grazing  incidence.  Again  the  AC  values  show  marked  asymmetry, 
similar  to  that  obtained  for  vertical  polarisation  at  6  GHz  with  the  0 
to  90*  measurements  up  to  11  dB  abovt  the  0  to  -90*  measurements.  The 
values  for  6  GHz  are  also  asymmetric  but  with  the  0  to  90*  vslues 
falling  some  10  dB  belov  the  0  to  -90*  measurements  (figure  15).  Once 
again,  theoretical  estimates  fall  well  below  measured  values. 

At  8  GHz  the  theoretical  and  experimental  vslues  agree  to  within  tZ  dB 
over  most  of  ths  range;  interestingly  the  exceptions  to  this  are  the 
minima  at  22*  and  48*  where  the  deeper  nulls  of  the  PO  solution  are  more 
suggestive  of  the  AC  data  than  the  GTD  and  UTD  solutions  which  have 
shallow  minima  some  B  dB  above  the  measured  vslues  (figure  16).  Near 
grazing  incidence  UTD  and  PO  estlmstss  are  well  below  AC  vslues  while 
GTD  values  are  well  above,  becoming  infinite  at  grazing  incidence. 

The  AC  data  at  10  GHz  show  significantly  greater  oscillatory  behaviour 
than  the  GTD  and  UTD  solutions,  though  the  Mpha»;"  of  the  variations  is 
in  rough  agrsemsnt  at  all  angles  (figure  17) .  Asymmetry  in  the  AC  data 
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suggests  sessurtMnts  are  affected  by  tha  problems  discussed  In 
Section  5.4, 

Measurements  at  12  GHz  are  In  slightly  better  agreement  with  the  GTD  end 
UTD  solutions  than  the  corresponding  vertical  polarise-  Ion  values  at 
this  frequency  but  variations  of  up  to  5  dB  occur  (figure  18).  The 
asymmetry  scatter  In  AC  data  again  calls  In  question  their  reliability. 

Summarising  the  results,  one  can  observe  that  at  the  lover  frequencies 
systematic  AC  data  structure  departs  significantly  from  theoretical 
predictions,  At  ths  Intermediate  frequencies,  8  end  10  GHz,  agreement 
Is  generally  very  good  particularly  for  vertical  polarisation.  Tha 
12  GHz  reaulta  are  laaa  reliable  for  the  reasons  mantioned  In 
Section  5.4. 

Except  at  the  lowest  frequency.  It  la  true  that  tha  PO  solutions  are 
satisfactory  for  representing  the  main  lobe  of  the  pattern,  near  normal 
Incidence;  at  higher  frequencies  reasonable  agreement  (♦  several 
decibels)  holds  for  the  second  and  third  lobes  as  well.  Of  course,  near 
the  minima  PO  must  not  be  used.  At  higher  frequencies  It  locates  the 
positions  of  the  minims  ressonahly  well  but  at  low  frequencies  it  falls 
even  to  do  that. 

The  GTD  solutions  have  a  singularity  at  gracing  Incidence,  a  failure  of 
the  method.  Apart  from  angles  near  edge  on,  the  GTD  and  UTD  solutiona 
are  In  better  agreement  with  AC  measurements  virtually  everywhere, 
though  at  low  frequanciea  discrepancies  of  6  to  16  dB  occur.  At  8  and 
10  GHz  the  GTD  and  UTD  solutions  ara  in  excellent  agreement  with  AC  data 
but  data  unreliability  at  12  GHz  precludes  any  extrusion  of  the  observed 
trend  to  that  frequency.  For  angles  near  grazing  Incidence  GTD  falls 
whereas  UTD  solutions  remain  finite;  agreement  between  UTD  solutions 
and  AC  data  is  reasonable  for  vertical  polarisation  but  not  so  for 
horizontal  polarisation,  UTD  solutions  being  much  lower  than  AC  data 
there. 

6.2  Rectangular  cylinder 

6.2.1  Vertical  polarisation 

For  the  10  cm  cylinder  the  AC  measurements  are  up  to  12  dB  greater  chan 
PO,  CTO  and  GTS  predictions,  though  s  merited  systematic  asymmetry  in  the 
experlmentsl  values  accounts  for  some  of  this  (figure  19).  The  expected 
ayaMtrv  about  45*  is,  however,  only  present  In  the  0  to  90*  sueasure- 
ments  wnich  are  the  ones  which  differ  most  from  ths  theoretical  curves. 
The  30  cm  cyllndsr  results  show  agreement  to  within  2  dB  everywhere 
except  at  the  sharp  minima  near  22*  and  68*.  where  GTD  and  UTD  are  too 
conservative  while  PO  exaggerates  ths  minima  (figurs  20).  At  50  cm 
cylinder  lengths  even  the  minima  are  accurately  predicted  by  GTD  and 
UTD,  with  excellent  agreement  at  all  angles  (figure  21). 

The  6  GHz  measurements  for  the  10  cm  cylinder  lie  well  above  theoretical 
predictions  at  intermediate  angles  hut  agree  reasonably  wall  near  0*, 
45*  and  90*  (figurs  22).  It  is  interesting  to  note  that  tha  AC  data 
show  slight  maxima  around  45*  as  predicted  by  PO  but  not  by  GTD  and  UTD. 
Ths  30  cm  cylinder  measurements  are  In  much  cloaar  agreement  with  ths 
GTD  and  UTD  calculations,  differing  by  at  most  3  dB  (figure  23).  PO 
consistently  underestimates  tha  ICS,  especially  near  minima.  For  the 
50  cm  cylinder  the  GTD  and  UTD  solutions  match  the  AC  data  equally  well, 
to  within  leas  than  the  difference  between  the  0  to  90*  and  0  to  -90* 
measurements  (figure  24).  The  PO  disagrees  as  before. 
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At  8  CHx  the  theoretical  results  for  the  10  cm  cylinder  match  well  out 
to  20*  where  they  all  underestimate  the  RCS  (figure  25).  Near  $  -  45* 
the  GTD  and  UTD  estimates  overestimate  and  predict  a  shallow  minimum; 
FO  Is  in  better  agreement  over  this  region  but  diecrepencles  of  up  to 
7  dB  are  present.  The  reverse  is  true  for  the  30  cm  cylinder 
(figure  26).  Near  $  *  45*  the  theoretical  results  are  too  low  but  the 
GTD  and  UTD  values  are  only  a  few  decibels  In  error.  The  nulls  at  11* 
and  22*  and  their  complementary  anglea  are  accurately  modelled  by  GTD 
and  UTD  while  PO  locates  them  correctly  but  overestimates  their  depths. 
Fox  the  50  urn  cylinder  GTD  and  UTD  provide  accurate  aolutlona  at  all 
tiapects  (figure  27).  A  slight  asyveatry  In  the  experimental  results 
suggests  an  equipment  effect. 

Turning  to  the  figures  for  10  GHx,  the  10  cm  cylinder  AC  measurements 
show  deep  nulls  at  37*  and  53*  (45*  1  8*)  which  are  not  adequately 
predicted  by  PO  and  are  only  hinted  at  by  GTD  and  UTD  (figure  28). 
Elsewhere  the  agreement  between  theory  and  experiment  Is  less  satisfac¬ 
tory  than  at  8  GHz.  The  30  cm  cylinder  results  are  in  much  better 
accord,  however..  with  GTD  and  UTD  estimates  within  2  dB  of  measured 
values  (figure  2$).  Apart  from  the  nulls,  the  PO  estimates  are  also 
good  everywhere  except  near  6  -  45*  where  they  depart  from  AC  data  by 
several  decibels.  Similar  remarks  can  be  made  about  the  50  cm  cylinder 
RCS  pattern  at  10  GHz  though  here  one  cennot  observe  any  superior 
performance  of  GTD  and  UTD  over  PO  near  p  “  45*  (figure  28). 

It  la  apparent  that  measured  values  of  tha  RCS  of  the  10  cm  rectangular 
cylinder  at  12  GHz  are  scattered  more  widely  about  the  theoretical 
curves  than  at  lower  frequencies  (figure  31).  The  unreliability  of  the 
12  GHz  AC  data  has  bean  noted  In  Section  5.4. 

At  a  cylinder  length  of  30  cm  tha  agreement  between  theory  and  experi¬ 
ment  is  certainly  not  aa  good  as  at  10  GHz.  Tha  shape  of  the  patterns 
agree  everywhere  but  deviations  of  3  or  4  dB  occur  for  4>  >  15*;  the 
scatter  In  AC  data  may  explain  this  excursion  (figure  32). 

The  50  cm  cylinder  results  are  similarly  slightly  less  Impressive  than 
at  lover  frequenclce  but  the  greater  asymmetry  in  the  anecholc  chimber 
data  may  account  for  this  (figure  33).  The  PO  solutions  in  this  case 
tend  to  underestimate  the  widths  of  the  lobes  for  $  >  20*. 

6.2.2  Horizontal  polarisation 

Tha  anecholc  chamber  measurements  of  the  10  cm  cylinder  at  4  GHz  ..re 
generally  In  batter  accord  vltb  GTD  and  DTD  solutions  thsn  was  the  case 
for  vertical  polarisation.  Tha  nulls  In  tha  pattern  are  much  mo.-e 
pronounced  for  horizontal  polarisation  (figure  34).  Thera  la  reasonable 
agreement  between  the  UTD  eolutlon  and  AC  data  in  tha  general  shape  o.: 
the  curve,  though  dlscrepenclea  of  several  decibels  are  observed  anc 
UTD' a  locations  of  the  nulls  are  out  by  ^<3*  from  AC  results.  The  same 
la  true  of  GTD  solutions  excapt  naar  6  “  0*  and  6  *  90*  where  GTD 
solutions  are  Infinite.  PC' s  prediction  of  the  main  lobes  agree*  well 
with  AC  data  but  it  falls  to  dcscriba  tha  behaviour  outside  the  main 
lobaa;  In  particular.  It  miases  the  lobe  at  45*  entirely. 

The  AC  date  for  the  30  cm  cylinder  are  marked  by  four  outliers  in  the  0 
to  90*  measurement e  (figure  35).  These  suggest  nulls  of  25  dB  but  tha  0 
to  -90°  and  theoretical  curves  era  in  agreement  that  tha  nulls  ars  only 
^15  dB.  Apart  from  these  outlier  points  tha  AC  data  la  In  fair  agree¬ 
ment  with  UTD  and  with  GTD  except  near  tha  singularities  at  0*  and  90*. 
Tha  FO  curve  la  within  1  dB  over  tha  main  lobaa  but  again  FO  falls  to 
predict  tha  lobe  at  45*. 
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At  *  cylinder  length  of  SO  on,  PO  again  represente  the  aein  lobes  to 
vlthln  1  or  2  dB  and  misses  the  lob*  at  AS*  (figure  36).  UTD  end.  apart 
fro*  the  effects  of  the  singularities  at  0*  and  90* >  GTD  are  In  good 
agreeaent  with  AC  data  on  the  shape  of  the  curve  but  AC  data  consis¬ 
tently  indicates  deeper  nulls  than  predicted  by  UTD  or  GTD. 

The  6  Gtis  Measurements  of  the  10  cm  cylinder  reveal  an  alaoat  global 
underestimation  of  RCS  by  the  theories,  ranging  from  a  few  decibels 
around  the  main  lobe  for  UTD  and  PO  to  10  dB  for  UTD,  10  dB  for  GTD  and 
15  dB  for  PO  near  45*  (figure  37),  GTD  again  predicting  infinite  fields 
in  the  main  lobes  at  0*  and  90*.  The  theories  fall  to  predict  curves 
that  agree  in  shape  with  that  of  measurements  between  30*  and  60* .  This 
is  uncharacteristic  of  the  GTD  and  UTD  solutions  thus  far.  While  some 
asymmetry  between  the  0  to  90*  and  0  to  -90*  measurements  is  evident  for 
$  >  75*.  the  theories  are  clearly  Inadequate  for  this  case.  It  is  worth 
noting  that  the  vertical  polarisation  results  showed  similar  underestl- 
Mtes  of  the  RCS  at  Intermediate  angles. 

The  30  cm  cylinder  data  are,  in  contrast,  quite  well  matched  by  UTD  and, 
apart  from  the  loba  at  45*,  by  PO  (figure  38).  GTD  besides  having 
singularities  at  0*  and  90*.  estimates  broader  Min  lobes  than  UTD,  PO 
and  AC  data. 

The  results  ror.  the  50  cm  cylinder  show  excellent  agreement  between  UTD 
and  the  AC  meuaurements  (figure  39).  Again  GTD  estiMtes  oroader  main 
lobes  thau  UTD,  PO  and  AC  data. 

At  8  GHz  the  theories  exaggerate  the  second  and  third  nulls  somewhat  but 
generally  the  10  cm  cylinder  RCS  pattern  is  quite  well  modelled 
(figure  40),  especially  compared  with  the  performance  of  the  theories  at 
4  and  6  GHz.  UTD  and  GTD  disagree  with  PO  In  the  location  of  the  second 
nulls  and  there  is  some  scatter  in  the  measurement*  making  location  of 
them  by  the  AC  data  uncertain.  UTD' a  estimates  of  the  first  nulls  agree 
very  wall  with  AC  data,  while  GTD  underestimates  them  somewhat  and  PO 
grossly  exaggerates  then  ss  it  does  the  second  nulls. 

The  GTD  and  UTD  calculations  for  the  30  cm  cylinder  predict  a  distortion 
of  the  second  lobe  near  22*  and  68*  which  does  not  appear  to  be  present 
In  the  AC  data  (figure  41).  Again  UTD  has  the  best  estlmatss  of  the 
first  nulla  while  both  UTD  and  GTD  underestlMte  the  second  nulla  and, 
as  before,  GTD  is  singular  at  0*  and  90*.  PO  exaggerates  the  first  and 
second  nulls  but  undersstlMtes  ths  lobes  at  30*  and  60*. 

The  situation  la  much  the  mm  for  the  50  cm  cylinder  (figure  42).  Ths 
distortion  of  ths  second  lobe  Is  again  apparent  In  the  GTD  and  UTD 
eolutlona  and  toms  suggestion  of  a  genuina  asymmetry  near  tba  paak  of 
the  lobe  la  present  In  the  AC  date.  The  Mtch  between  theory  and 
cquipmant  is  quits  good  slsewhara,  within  2  to  3  dB  evaryvhere  for  UTD 
excapt  for  the  second  null,  within  2  to  3  dB  for  GTD  except  for  the 
first  end  second  nulls  and  near  0*  and  90*  and  within  2  to  4  dB  for  PO 
except  et  the  nulls. 

The  anecholc  chamber  Muureienti  et  10  GHz  for  the  10  cm  cylinder  are 
somewhat  erratically  distributed  with  gross  variations  bstveen  tha  0  to 
90*  and  0  to  -90*  pattarn  -  in  toms  cases  they  differ  by  18  dB 
(figure  43).  Characterisation  of  the  Min  lobes  (sway  from  0*  end  90* 
for  GTD)  la  good  but  the  scatter  of  the  measurement  Mkes  further 
interpretation  difficult. 
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More  sense  can  be  made  of  the  raault  for  th«  30  cm  cylinder  (figure  44). 
Here  Che  agreeaent  between  theory  and  experiment  le  obviously  very  good, 
typically  11  to  12  dB,  with  only  Isolated  exceptions.  For  Instance,  P0 
appears  to  aialocate  the  null  near  29*,  placing  It  at  27*  whereas  GTD 
and  UTD  matches  the  measured  data.  GTD  underestimates  the  first  nulls 
while  UTD  overestlaates  the  second  nulls. 

The  30  en  cylinder  likewise  shows  very  good  agreement  between  UTD  and 
aeasureaents  though  there  is  sore  scatter  in  the  measurements  for 
$  >  60*  (figure  45).  Away  froa  0*  and  90*,  GTD,  apart  from  under¬ 
estimating  the  first  nulls,  is  equally  Impressive.  PO  again  fails  to 
predict  the  correct  behaviour  for  20*  to  30*  and  its  Image  for  70*  to 
60*.  Apart  from  this,  it  too  agrees  well  with  measurements. 

The  measurements  at  12  GHz  for  the  10  cm  cylinder  (figure  46)  show  a 
fair  amount  of  scatter  between  20*  to  70*  as  la  the  case  for  the  10  GHz 
measurements  of  the  10  cm  cylinder.  Outside  of  the  20*  to  70*  region, 
agreement  between  theories  and  aeasureaents  is  good  (away  from  0*  and 
90*  for  GTD).  In  the  20*  to  70*  region  PO ' s  prediction  departs  signifi¬ 
cantly  from  that  of  GTD  and  UTD;  la  particular  the  peaks  predicted  by 
GTD  and  UTD  at  30*  and  60*  ara  missed  by  PO  which  placee  deep  nulla 
there  Instead.  Deaplta  tha  poor  measurements  here,  more  confidence  is 
to  be  placed  in  the  GTD  and  UTD  results. 

The  results  for  the  30  cm  cylinder  show  vsry  good  agreement  between  UTD 
and  AC  data  apart  from  an  overastlmation  of  tha  peaks  at  30*  and  60*  by 
UTD  (figure  47).  The  same  la  true  of  GTD  away  from  0*  and  90*.  PO 
overestimates  the  flrct  two  nulls  from  the  mein  lobes  and  rsplacas  the 
peaks  at  30*  and  60*  by  daep  nulls  as  before;  apart  from  that.  It 
agraas  wall  with  maaauraments . 

The  50  cm  cylinder  measurements  have  mors  scatter,  particularly  froa  70* 
to  90*  (figure  48) .  Good  agreeaent  la  obtained  between  UTD  and  AC  data 
and,  away  from  0*  and  90*,  between  GTD  and  AC  data,  especially  with  the 
0  to  90*  measurements.  P0  suffers  from  the  same  deficiencies  as  men¬ 
tioned  for  the  30  cm  cylinder.  Nevertheless,  the  height  of  the  peaks  as 
predicted  by  PO  matches  that  of  measurements  better  than  that  of  GTD  and 
UTD  which  overestimate  them  particularly  away  from  tha  main  lobes 
directions. 

6.3  Comparison  with  momant  method  solution 

The  results  obtained  In  this  paper  for  the  flat  plats  and  finite 
rectangular  cylinder  using  the  normalisation  schema  through  aquation  (10) 
must  Involve  termination  errors,  the  significance  of  which  must  Increase 
with  decreasing  plate  size  and  cylinder  length.  Comparison-  with 
experimental  results  attaats  to  this  expected  trend  in  the  cylinder  case, 
agreement  Is  good  down  to  a  length  of  21  but  In  the  case  of  the  flat  plata 
it  la  hard  to  aaaess  the  accuracy  of  tha  solutions  obtained  ee  there  Is 
much  more  scatter  In  the  AC  data  ss  has  been  mentioned  In  Section  5.4. 

Xn  view  of  the  unreliability  of  tha  AC  data  at  low  power  return,  the  PO  and 
UTD  results  are  compared  here  with  moment  method  (KM)  solutions  of  the 

2b  2a 

Integral  equation  for  the  flat  plata  at  4  GHz,  that  is,  for  ~  "  ~  •  1.3. 

Comparison  is  not  made  with  GTD  as  it  gives  tha  asms  results  as  UTD  away 
from  grazing  incidence.  Tha  KM  solutions  were  obtained  using  tha  NIC  wire 
grid  modelling  code  available  at  DXCS.  Agreement  between  PO  and  KM 
solutions  la  only  good  for  the  main  lobes,  that  Is,  user  tha  apacular 
directions  (figures  49  and  50).  0s  tha  othmr  band  agreement  between  UTD 
end  !M  eolations  at  4  GHz  is  very  good  -  within  l  dB  up  to  25*  and  within 
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2  dB  bayond  25*  for  vtrtical  polarisation  (figure  49),  within  1  dB  up  to 
40*  and  within  4  dB  bayond  40*  for  horizontal  polarisation  (flgura  50). 
Thaaa  raaulta  indicate  that  tarnlnatlon  affacta  are  not  significant  even 
for  plata  lengths  as  small  as  1.3X. 

Figures  51  and  52  compare  DTD,  PO  and  KM  solutions  for  the  flat  plata  at 
2b  2a 

2  GHz,  that  is,  for  “  “  “  “  0.7,  for  vertical  and  horizontal  polarisation 

respectively.  At  this  frequency,  aven  PO' a  prediction  of  the  specular 
return  from  the  plate  la  several  decibels  out  froa  that  of  MM.  DTD's 
solution  deviates  froa  that  of  M(  by  as  such  as  2  dB  up  to  40*  and  8  dB 
beyond  40* ,  agreeaent  being  slightly  better  for  horizontal  polarisation. 
Thus  if  a  Halt  has  to  be  placed  on  the  aethod  of  DTD  and  CTD  (away  froa 
its  singularities)  and  equation  (10)  in  the  prediction  of  the  RCS  of  the 
flat  plate  it  la  for  a  plata  of  aize  MX . 

The  UTD  solution  for  the  rectangular  cylinder  has  been  shown  to  be  good 
down  to  cylinder  width  “  1.3X  and  length  2. OX.  Coaparison  of  results  with 
MM  solutions  would  reveal  the  lower  Halt  of  validity,  independent  of  data 
reliability.  This  coaparlaon  la  scarcely  feasible  as  can  be  appreciated  by 
noting  that  the  MM  solution  for  the  1.3X  flat  plate  took  30  ain  of  CPU 
tine.  In  view  o£  this  no  atteapt  was  made  to  obtain  MM  solutions  tor  the 
rectangular  cylinder. 

6.4  Deficiencies  of  the  theories 

5.4.1  Physical  Optics 

It  is  a  common  feature  of  the  PO  solutions  that  they  poaaess  deep  nulls 
in  RCS  which  are  not  realised  in  the  experlaental  aeasureaents .  This 
effect  could  have  been  eliminated  froa  the  figures  by,  for  example, 
convolving  the  PO  solutions  with  a  windowing  function.  More  conve¬ 
niently,  smoothing  could  have  been  done  by  weighting  in  the  transform 
domain.  No  such  processing  has  been  perforated  here  to  emphasise  the 
fact  that  simple  methods  such  as  PO  do  yield  unpbyslcal  results. 

In  general,  PO  proves  to  be  accurate  in  the  prediction  of  the  main  lobes 
of  the  RCS  patterns,  accuracy  falls  away  for  aspect  angles  away  froa 
theae  regions.  This  arises  froa  the  fact  that  the  surface  field  is 
taken  to  be  the  geoaetric  optles  field.  A  aethod  developed  by 
Uflatssv(ref . 10)  known  as  the  physical  theory  of  diffraction  corrects 
for  this  surface  field  in  PO  by  including  an  additional  term  in  the 
induced  surface  current  -  this  '  nonvalform'  component  is  due  to  pertur¬ 
bations  created  by  the  departure  of  the  surface  froa  an  infinite  plane 
and  may  arise  froa  the  curvature  of  the  surface  at  a  shadow  boundary  or 
froa  geometrical  discontinuities  such  as  edgaa. 

The  backscattered  RCS  as  predicted  by  PO  is  polarisation  independent,  a 
result  of  the  tangent  plana  approximation.  This  is  a  fundamental 
weakness  of  PO  in  general  though  it  can  yield  polarisation  dependent 
results  in  the  case  of  ecatterers  with  multiple  reflections,  eg  ref¬ 
erence  3. 

6.4.2  Geoaatrlcal  theory  of  diffraction  and  uniform  gaoaatrlcal  theory 
of  diffraction 

CTD  falls  in  the  transition  regions  adjacent  to  shadow  and  rsflaction 
boundaries  where  the  solution  becomes  Infinite  aod  discontinuous .  This 
limitation  is  overcome  by  the  uniform  GID(UTD)  developed  by  Kouyoumjlan 
and  Pathek(ref *8)  and  by  the  uniform  field  introduced  by 
Abluwalla  at  sl(ref . 11 , 12) .  Both  approaches  yield  high-frequency  fields 


011.-0344-1* 


34  - 


that  are  continuous  at  shadow  and  reflection  boundaries  and  that  reduce 
to  the  field*  predicted  by  GTD  away  from  euch  transition  regions. 


The  singularities  that  arise  In  the  GTD  solution  for  the  strip  and  the 
cylinder  appear  in  the  aultlple  diffraction  terms  at  shadow  boundaries. 
They  arc  a  consequence  of  using  Keller's  plane  wave  diffraction  coeffi¬ 
cients  In  ths  treatment  of  higher  order  diffraction.  When  the  higher 
order  diffraction  terms  are  formulated  In  terms  of  cylindrical  wave 
diffraction,  as  In  UTD,  the  results  are  finite. 


It  la  Interesting  to  note  that  the  singly  diffrected  GTD  or  UTD  field 
for  the  strip  (equation  (16)  for  vertical  polarisation  and  equation  (39) 
for  horltontal  polarisation).  Is  finite  at  4  ■  0*  despite  the  fact  that 
the  individual  terms  are  infinite  there.  The  same  Is  observed  to  hold 
for  the  rectangular  cylinder  (equation  (53)  for  vertical  polarisation 
and  aquation  (59)  for  horltontal  polarisation),  at  p  ■  0*  and  $  -  90*. 
This  Is  a  phsnomsnon  that  occurs  frsqusntly  In  CTD/UTD  snslysls  of 
rsdlstlon  snd  scattering  problems.  It  wss  noted  by  .Temee(ref .6)  that 
the  endpoint  Infinities  would  cancel  to  yield  s  finite  result  If  the 
aperture  had  a  symmetrical  field  distribution,  though  in  general  this 
reault  is  not  expected  to  be  the  same  as  that  obtained  by  e  direct 
evaluation  of  the  aperture  field  Integral. 

A  fundamental  limitation  of  GTD  and  UTD,  In  common  with  other  rey 
techniques,  is  the  prediction  of  Infinite  fields  at  caustics.  This  doss 
not  correspond  to  physical  reality  snd  caustic  corrections  have  to  be 
Introduced.  No  general  techniques  valid  for  an  arbitrary  body  are 
available  for  finding  caustic  corrections;  this  Involves  generating 
matching  functions  thst  provide  a  smooth  transition  Into  a  finite  value 
at  ths  c*ustlcs(ref.l3,14,i5).  Tbs  GTD  and  UTD  edge  diffraction 
concepts  are  suitable  for  ths  analysis  of  geoatetries  where  diffraction 
appears  to  coma  from  a  single  point  or  group  of  single  Isolated  points 
along  an  edge.  At  $  -  0*  for  the  flat  plate,  and  #  •  0*  and  90*  for  the 
finite  rectangular  cylinder,  there  la  e  congruence  of  backscattered 
raye,  that  la,  caustics  occur.  Thus  GTD  and  UTD  fall  to  predict  the 
field  In  the  backscattered  directions.  In  these  cease  it  la  convenient 
to  use  the  physical  optics  approximation  to  correct  the  scattered  field 
determined  by  GTD  end  UTD;  this  caustic  correction  technique  la  ueeful 
for  the  flat  plate  structures  treated  here.  Another  technique  for 
caustic  correction  Is  to  employ  equivalent  electric  and  magnetic  edge 
currants  derived  from  the  GTD /UTD  diffracted  field  away  from  the 
caustics.  These  equivalent  currents  are  substituted  in  the  radiation 
Integral  to  give  the  caustic  field. 


In  this  psper  ths  fist  piste  end  finite  rectangular  cylinder  have  not 
been  treated  aa  three  dimensional  scattering  targets;  rather  the 
results  were  obtained  by  using  a  normalisation  scheme  (equation  (10))  to 
reduce  the  two  dimensional  Infinite  strip  and  infinite  cylinder  solu¬ 
tions  to  onas  appiicabla  to  ths  finite  problems.  Interestingly,  ths 
solutions  thus  obtained  in  the  backscattered  directions  corresponding  to 
caustics  for  the  finite  target  problem,  yield  finite  values  Chat  agra* 
with  PO  even  though  ths  geometrical  optics  reflection  terms  were  not 
included.  This  effect  has  been  observed  by  Jaaes(ref.6)  In  his 
treatment  of  radiation  from  s  parallel-plat#  waveguide  as  has  the 
cancellation  of  the  singularities  In  the  singly  diffracted  field  terms 
mentioned  earlier. 


- - 
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7.  CONCLUSIONS 

PO  proves  to  be  very  accurate  in  the  neighbourhood  of  specular  reflections, 
that  is,  in  the  region  of  the  main  lobes  of  the  RCS  curves,  at  leant  down  to 
plate  sizes  2a  'V  1.3X.  Its  perforaance  deteriorates  with  aspect  angle  away 
from  the  specular  direction  where  polarisation  dependent  effects  become 
evident  but  are  not  predicted  by  PO,  though  the  range  of  reliable  prediction 
increases  with  plate  size.  It  should  be  noted,  hov^ver,  that  apart  from  a 
substantial  overestimatior.  of  the  depths  of  the  nulls,  PO's  domain  of  applic- 

2a 

ability  is  extended  to  the  next  sidelobes  along  with  increasing  -p 

Singularities  exist  in  the  GTD  solution  for  both  polarisations  at  grazing 
Incidence  on  the  flat  plate.  Away  from  the  neighbourhood  of  these  singular¬ 
ities  ,  the  GTD  solution  for  the  flat  plate  agrees  with  the  UTD  solution.  GTD 

2a 

and  UTD  results  for  the  flat  plate  improve  as  —r  increases.  The  AC  data  at 

2a  A 

the  smaller  values  of  are  erratic  but  comparison  with  MM  solutions  has 

demonstrated  the  accuracy  of  UTD  solution  down  to  plate  dimensions  of  one 
wavelength. 

An  ideal  thin  flat  plate  should  yield  a  zero  RCS  at  grazing  incidence  for 
horizontal  (that  is,  transverse)  polarisation  in  th«  limit  of  zero  thickness, 
yet  the  measured  RCS  for  horizontal  polarisation  does  not  diminish  to  zero  at 
grazing  incidence  as  expected.  This  suggests  that  the  plate  is  not  a  "thin" 
plate  as  assumed  or  that  an  unexpectedly  significant  contribution  comes  from 
vertex  diffraction.  The  thickness  or  the  edge  ranges  from  0.01X  at  4  GHz  to 
0.04X  at  12  GHz.  do  exact  solution  is  available  for  diffraction  from  an  edge 
with  finite  thickness,  however  references  16  and  17  present  solutions  for  a 
thick  edge  modelled  as  two  90”  wedges.  The  application  of  this  model  to  the 
present  case  has  not  been  attempted. 

The  GTD  solution  for  the  rectangular  cylinder  for  horizontal  polarisation  has 
singularities  at  -  C  and  90*.  Away  from  the  neighbourhood  of  these 
singularities  close  agreement  is  obtained  between  GTD  and  UTD  solutions.  The 

2a 

results  for  the  rectangular  cylinder  at  4  GHz  work  well  down  to  1.3  for 

the  30  cm  and  50  cm  cylinders  where  termination  errors  are  not  significant. 
Agreement  of  the  UTD  solution  with  measured  results  is  definitely  not  good  for 
the  10  cm  cylinder,  particularly  at  the  lower  frequencies  of  4  and  6  GHz; 
data  reliability  is  questionable  in  these  cases  and  in  the  absence  of  compari¬ 
son  wirh  accurate  numerical  solutions  of  the  exact  integral  equations  for 
scattering  from  the  cylinder,  it  is  difficult  to  assess  the  performance  of  the 
method  for  these  cases.  It  could  well  be  that  termination  errors  arising  from 
the  assumptions  implicit  in  equation  (10)  are  more  significant  for  the 
cylinder  of  length  1.3X  than  for  the  flat  plate  of  length  1.3X.  It  seems  more 
likely,  however,  that  the  poor  agreement  for  the  10  cm  cylinder  is  due  to 
errors  in  the  anechoic  chamber  measurements,  as  in  the  case  of  the  flat  plate, 
and  that  UTD  with  equation  (10)  is  a  good  description  of  scattering  from  a 
finite  cylinder  down  to  cylinder  lengths  of  the  order  of  a  wavelength. 
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Figure  IS.  RCS  for  flat  plate  at  &  GHz  (horizontal  polarisation) 


or  flat  plate  at  8  GHz.  (horizontal  polarisation) 
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Figure  18.  RCS  for  flat  plate  at  12  GHz  (horizontal  polarisation) 
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Figure  19 


Figure  19.  RCS  for  10  cm  cylinder  at  4  GHz  (vertical  polarisation) 
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for  30  cm  cylinder  at  4  GHi-  (vertical  polarisation) 
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Figure  22.  RCS  for  10  cm  cylinder  at  6  GHi  (vertical  poiarisat 
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Figure  24 
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Figure  25 


Figure  2S.  RCS  for  10  cm  cylinder  at  6  GHz  (vertical  polarisation) 
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Figure  26.  RCS  for  30  cm  cylinder  at  8  GHz  (vertical  polarisation) 
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Figure  27 
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:s  for  50  cm  cylinder  at  10  GHz-  (vertical  polarisation) 
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Figure  31 


figure  32.  RCS  for  30  era  cylinder  at  12  GHz  (vertical  polarisation) 
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Figure  33 
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Figure  33.  RCS  for  50  cm  cylinder  at  12  CHz  (vertical  polarisation) 


figure  34.  ItCS  for  10  cm  cylinder  at  4  GHz  (horizontal  polarisation) 
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Figure  35 
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50  cm  cylinder  at  A  GHz  (horizontal  polarisat 


figure  38.  RGS  for  30  cm  cylinder  at  6  GHz  (horizontal  polarisation) 
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Figure  40.  HCS  for  10  era  cylinder  at  8  Cllz  (horizontal  polarisation) 
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Figure  42 


igurc  44.  RCS  for  30  cm  cylinder  at  10  GHz  (horizontal 
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16  |  SUMMARY  OR  ABSTRACT: 

I  (if  this  is  security  classified,  the  announcement  of  this  report  will  be  similarly  classified) 


The  method  of  Physical  Optics,  the  Geometrical  Theory  of  Diffraction 
and  the  Uniform  Geometrical  Theory  of  Diffraction  are  applied  here  to 
electromagnetic  scattering  from  flat  plates  and  rectangular  cylinders. 
The  results  are  compared  with  experimental  measurements  to  establish 
the  domains  of  validity  of  the  theoretical  techniques.  * 
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